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Abstract. Wc investigate the semigroup structure of bosonic 
Gaussian quantum channels. Particular focus lies on the sets of 
channels which are divisible, idempotent or Markovian (in the sense 
of cither belonging to one-parameter semigroups or being infinites- 
imal divisible). We show that the non-compactness of the set of 
Gaussian channels allows for remarkable differences when compar- 
ing the semigroup structure with that of finite dimensional quan- 
tum channels. For instance, every irreversible Gaussian channel is 
shown to be divisible in spite of the existence of Gaussian channels 
which are not infinitesimal divisible. A simpler and known con- 
sequence of non-compactness is the lack of generators for certain 
reversible channels. Along the way we provide new representations 
for classes of Gaussian channels: as matrix semigroup, complex 
valued positive matrices or in terms of a simple form describing 
almost all one-parameter semigroups. 



1. Introduction 

A quantum channel describes the input-output relation of a quan- 
tum mechanical operation. Mathematically, it is described by a com- 
pletely positive map which is trace-preserving (in the Schrodinger pic- 
ture) or identity-preserving (in the Heisenberg picture). Considering 
equal input and output spaces we can concatenate quantum channels 
with themselves or other channels. Clearly, such a concatenation is 
again a valid quantum operation, so that the set of quantum chan- 
nels forms a semigroup. Inverses only belong to this semigroup if the 
channels describe unitary evolution. Apart from this subgroup of re- 
versible channels various subsets can be distinguished according to their 
semigroup properties; for instance (i) channels which are elements of 
one-parameter semigroups, i.e., physically speaking solutions of time- 
independent Markovian master equations, (ii) channels which are infin- 
itesimal divisible, e.g., solutions of time-dependent Markovian master 
equations, (iii) channels which are divisible into others in a non-trivial 
way, (iv) channels which are indivisible, and (v) channels which are 
idempotent. 
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Whereas the characterization of one-parameter semigroups goes mainly 
back to the seventies [H [2] the distinction of the above sets has been 
addressed more recently in [1] and methods to decide membership in (i) 
have been provided in [5]. The classical counterpart known as imbed- 
ding problem for Markov chains was exhaustively studied in probability 
theory (see [3] and the references therein.) 

The characterization of the above sets seems to become vastly more 
complex as the dimension d of the considered system is increased: 
whereas for qubits {d = 2) essentially everything is known, the problem 
of deciding membership in (i) turns out to be NP-hard with increasing 
d 1]. 

The present article is devoted to the study of the semigroup prop- 
erties of bosonic Gaussian channels (i.e., 'quasi-free' maps) — a class 
where the underlying Hilbert space is infinite dimensional. The re- 
striction to Gaussian channels is motivated by their physical relevance 
(they model, for instance, optical fibers and occur all along with qua- 
dratic interactions) and it is suggested by the mentioned complexity 
issue. Note that this restriction has two flavors: when we ask whether 
a channel is in one of the above sets we do not only restrict the chan- 
nel under consideration, we also restrict the involved one-parameter 
semigroups or factorizations to within the Gaussian setting. In this 
way we can escape from infinite dimensional Hilbert space into finite 
dimensional phase space and formulate everything in terms of finite 
dimensional matrix analysis. 

At this point one might expect that the basic picture of the set of fi- 
nite dimensional quantum channels carries over to the Gaussian setting. 
There are, however, crucial differences, for instance: whereas the set of 
quantum channels in finite dimensions is compact, the set of Gaussian 
channels (even though having a finite dimensional parametrization) is 
not; every irreversible Gaussian channel turns out to be divisible in 
spite of the existence of Gaussian channels which are not infinitesimal 
divisible etc. 

The article is organized as follows: Secj2] introduces the basic no- 
tation and SecE] shows that the set of Gaussian channels is indeed a 
matrix semigroup. In SecJH which is mainly provided for complete- 
ness, we review the reversible case and emphasize the fact that not 
every canonical transformation has a generating quadratic Hamilton- 
ian. Sec|5] deals with one-parameter semigroups, for which a simple 
representation is provided, and SecE] has a closer look at infinitesi- 
mal divisible channels. In SecJ7]we prove that all irreversible Gaussian 
channels are divisible by exploiting a simple mapping from the set of 
Gaussian channels to the cone of complex positive matrices. Finally, 
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in Secs J8|9] idempotent and gauge- covariant channels are investigated 
before SecHU] concludes with some open questions. 

2. Gaussian channels 

Let Qj,Pj, j = l,...,n, be the canonical operators satisfying the 
the canonical commutation relations {h = 1) 

[Qj , Pfc] = i5jk , [Qj,Qk] = [Pj , Pk] = . 

We will use a notation R = (Qi, -Pi, . . . , Qn, PnY ^"^^ ^ach ^ G M^", 
we define = e*^ Here we have denoted 

a = an = ®'i=i<yi , en = 

The unitary operators are called Weyl operators and they corre- 
spond to displacements in phase space. 

A Gaussian channel is a quantum channel which maps Gaussian 
states into Gaussian states [711H1I9]- The mathematical form of Gaussian 
channels is best described in the Heisenberg picture when we look at 
their action on the Weyl operators. A Gaussian channel then corre- 
sponds to a mapping 

where X, Y are real 2n x 2?7,-matrices. Here and thereafter without loss 
of generality we restrict to channels which map zero-mean states into 
zero-mean states. Complete positivity (cp) imposes a constraint on the 
matrices X, Y, which can be written as 

Y >i{a- XaX^) . (1) 

From now on, we will identify the Gaussian channels with the pairs 
{X, Y) of real 2n x 2n-matrices satisfying ([T]). We denote by Q the set 
of all Gaussian channels. 

It is often useful to depict a Gaussian channel by its action on the 
first and second moments of a quantum state. We denote the first 
moments by a vector d G M^" whose components are expectation values 
dk = (Rk) and we define the covariance matrix as 

= {{Rk — dk, Ri — di}+) . 
A Gaussian channel (X, Y) then acts as 

d 1—^ Xd , 

r h-> XTX^ + Y. 
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Example 1. The preparation of a Gaussian state is a simple instance 
of a Gaussian channel. A Gaussian channel with X = has an input- 
independent output state with covariance matrix Y. The cp-condition 
then reduces to Y + ia > 0, which is nothing but the condition for 
y to be a valid covariance matrix. 

Generally speaking, the Y contribution in a Gaussian channel (X, Y) 
can be regarded as noise term. It follows from ([T]) that y > 0. If y = 0, 
then the cp-condition ([T]) implies that XaX^ = a, meaning that X is 
an element of the real symplectic group Sp{2n). The group of unitary 
Gaussian channels is therefore identified with Sp{2n). Let us notice, 
however, that generally X can be any real matrix, as long as sufficient 
noise is added (i.e., Y is large enough). 



Concatenation of two Gaussian channels is again a Gaussian channel. 
Hence, the set of Gaussian channels forms a semigroup. The semigroup 
product is given by 



and the identity element is (1,0). 

Let us recall that we are using the Heisenberg picture and in the 
Schrodinger picture the order of the product is opposite to the order of 
application of channels. For instance, if (Xi,Yi) and (^2,12) describe 
optical fibers, then the product (Xi,Yi) ■ (X2,F2) corresponds to the 
channel in which the signal first goes through (X2, Y2) and then through 



The semigroup product can also be written as an ordinary matrix 
product. For each X,Y E M2„(M), we denote x = X (S) X and y G 
M^" (g) is the column vector defined through the condition 



3. Semigroup product 



(Xi, Fi) ■ (X2, Y2) = (X1X2, Fi + X^Y2Xf) 



(2) 



{X,,Y,). 



(ei^ejly) = {ei\Yej) . 
Then the mapping tt : Q ^ M4n2+2n+i(IR) with 



(3) 




is an injective homomorphism. We conclude that the semigroup of 
Gaussian channels is a matrix semigroup — a useful property when it 
comes to the discussion of one-parameter semigroups. 
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4. Reversible channels 

A channel is called reversible if it has an inverse which is also a chan- 
nel (i.e., we mean physically reversible as opposed to mathematically 
invertible). Reversible channels are exactly the unitary channels, and 
a Gaussian channel {X, Y) is reversible iff X e Sp{2n) and Y = 0. 

For completeness we briefly discuss the group structure of reversible 
channels, which already exhibits some interesting features. The sym- 
pletic group Sp{2n) is a non-compact connected Lie group whose Lie 
algebra sp(2n) is given by all real matrices s such that (sa)^ = sa. 
The exponential map 

exp : sp(2n) Sp{2n) , s (4) 

is not surjective - a common feature of non-compact Lie groups. 

In physical terms, the lack of surjectivity of the exponential map 
means that there are canonical transformations S e Sp{2n) for which 
there is no 'Hamiltonian matrix' s G sp(2n) generating them via S* = e*. 
Clearly, if we consider the corresponding unitary Us acting on Hilbert 
space (i.e., an element of the metaplectic representation of Sp{2n) [TO]), 
then there is always a Hamiltonian H such that Us = expiH. Such a 
Hamiltonian may be obtained from the spectral decomposition oi Us- 
The point is, however, that H cannot be a quadratic expression in the 
canonical operators if S' ^ exp(3p(2n)). 

A necessary condition for S G exp(sp(2n)) is that S has a real log- 
arithm. We recall the following standard result from matrix analysis 

m- 

Proposition 1. A real matrix X G M„(M) has a real logarithm L G 
M„(M) (i.e. X = e^) iff X is non-singular and the Jordan blocks of X 
corresponding to negative eigenvalues have even multiplicities. 

The characterization of the set of exp(5p(2?7.)) goes back to Williamson 
|12j . For our purposes, the following partial characterization should 
suffice. 

Proposition 2. Let S G Sp{2n). If — 1 is not an eigenvalue of S, then 
S G exp(sp(2r2)) iff S has real logarithm. 

A simple consequence of Prop. [2] is that if S* G Sp{2n) is positive, 
then S G exp(5p(2n)). 

An important subgroup K{2n) of Sp(2n) is formed by those matrices 
which are in addition orthogonal, i.e., K{2n) = Sp{2n) fl 5*0(277.). The 
corresponding maps are called passive transformations as they preserve 
the number of particles [10]. The subgroup K(2n) is a (maximal) com- 
pact subgroup of Sp{2n) and it is isomorphic to U{n). Consequently, 
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the exponential map from the Lie algebra ^(2n) to K{2n) is surjective. 
That is, every S G K{2n) has the property that S G exp(^(2ra)) C 
exp(sp(2n)). 

Using the Euler decomposition (cf. [10]) each S G Sp(2n) can be 
written as a product S = K1DK2, where Ki,K2 G K(2n) and D G 
Sp{2n) is a diagonal matrix of the form D = diag((ii, l/rfi, . . . , 1/ dn) 
with di,...,dn > 0. The diagonal matrix D describes single-mode 
squeezings. Thus, every reversible Gaussian channel is a concatenation 
of passive transformations and single-mode squeezings. 

From the previous discussion we arrive at the following conclusion. 

Proposition 3. Let S G Sp{2n). There are 81,82 G exp(sp(2?T,)) such 
that 8 = 8182- 

Proof. Since 8 has an Euler decomposition 

8 = K1DK2 = KiDKjK^^K2 . 

it is a product of a positive 81 := KiDKj and an orthogonal 82 = 
K^^ K2 symplectic matrix. Both 5*1, 5*2 G exp(sp(2n)), as discussed 
earlier. □ 



5. One-parameter Gaussian semigroups 

5.1. General form. By one-parameter semigroup of Gaussian chan- 
nels we mean a family of Gaussian channels, parametrized by M"'', which 
satisfies the following conditions: 

(i) continuity: Xt and Yt depend on t G M"*" in a continuous way 

(ii) semigroup property: XtXs = Xt+s and Yt+s = Yt -\- XtYgXf . 

(iii) connected to the identity: Xq = 1 and Yq = 

As we have seen in Section [3], the semigroup of Gaussian channels 
is a matrix semigroup. It follows that every one-parameter semigroup 
of Gaussian channels has a generator [13]. In particular, the mappings 
t ^ Xt and t ^ Yt are differentiable. The following characterization 
(although in a slightly different form) has been derived in [H]. 

Proposition 4. A family (Xj, Yt)t>o of Gaussian channels forms a 
one-parameter semigroup iff there exists real matrices A, B, H with iA+ 
B>0,A'^ = -A, = H such that 

Xt = e*(^-^)% (5) 

Yt = 2 [ Xja^BaXs ds . (6) 
Jo 



THE SEMIGROUP STRUCTURE OF GAUSSIAN CHANNELS 7 



In the Heisenberg picture this corresponds to an evolution of any 
observable O governed by a master equation dtO = C{0) with a Liou- 
villian 

c{o) = z[h,o] + j2l:ol-1{l:L,ou, (?) 

a 

= 2 ^ HkiRkRi , (8) 

kl 

La = ^La,kRk, with B + iA = L*L. (9) 

k 

A simple consequence of this characterization is the following. 

Corollary 1. Let X G M2„(M). The following conditions are equiva- 
lent: 

(i) There exists a matrix Y G M2„(M) such that {X,Y) is an ele- 
ment of a one-parameter semigroup of Gaussian channels. 

(ii) X is non-singular and the Jordan blocks of X corresponding to 
negative eigenvalues have even multiplicities. 

Proof. By Prop{T] the condition (ii) on X is nothing but the existence 
of a real logarithm L which is clearly necessary for X to occur in an 
element of a one-parameter semigroup of Gaussian channels. Hence, 
(i) implies (ii). 

In order to see the other direction, suppose that (ii) holds and let L be 
a real logarithm of X. Let us decompose La^ = A—H into a symmetric 
part H = and anti-symmetric part A = —A'^, respectively. Then 
there is always a 5 G M2„(M) (e.g. B = \\A\\^ 1) such that iA + B >0 
and we can construct a one-parameter semigroup of Gaussian channels 
by following the characterization in Prop Si □ 

Prop, m gives a complete but rather cumbersome characterization 
of one-parameter semigroups of Gaussian channels. In particular, the 
appearing integral might complicate further use of the characteriza- 
tion. Fortunately, almost all generators of such semigroups allow for a 
simpler representation discussed in the next subsection. 

5.2. Simple form. Suppose that {Xt}t>o is a semigroup and fix a real 
symmetric matrix y. Then by setting 

Yt = y- XtyxJ (10) 

the semigroup property for (X^, Ft)t>o is satisfied. The cp-condition 
now reads 

y - la > Xtiy - ia)Xj . (11) 
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Example 2. An amplification channel is of the form X = y/r] 1, 
Y = (?7 — 1) 1 for some t] G (l,oo). Amplification channels form a 
one-parameter semigroup. Namely, 

X, = e*l, Yt = {e^'- 1)1. 

This is of the simple form ( fTOl) with y = —1. 

Proposition 5. A one-parameter Gaussian semigroup {Xt, Yt)t>o is of 
the simple form (fTOj) if the operator {A — H)a in Prop. [7] (i.e., the 
generator of Xt) does not have a pair of eigenvalues of the form ±A. 

Proof. The general form ([6]) gives Yo = 2 a^Ba = B. On the other 
hand, the simple form (ITOl) leads to 

Yo = Ay + yA^, 

where A = {H — A)(t. Hence, in order for {Xt, Yt)t>Q to be of the form 
(JTO|l . we need to find y such that 

Ay + yA^ = B. (12) 

The linear equation (fT2|) can be written as 

(l®i^ + i® 1)13^) = \B) 

Hence, if 1 ® A^ + A ® 1 is invertible, then we have a solution for y. 

The eigenvalues of 1®A^ + A®1 are of the form Aj + Aj, where Aj, Xj 
are eigenvalues of A (this can be seen e.g. using Schur upper-triangular 
form for A and A^). Therefore, if A has the property that the sum of 
any two of its eigenvalues is nonzero, the invertibility of 1 AJ' + A ® 1 
follows. 

If 3^ is a solution then y'^ and hence \{y + 3^"^) are again solutions. 
Therefore, y can be chosen symmetric. □ 

PropjS] implies that almost all one-parameter semigroups admit a 
representation of the simple form in Eq. ffTOl) . Moreover, it shows that 
we can approximate any one-parameter Gaussian semigroup with a one- 
parameter semigroup of the simple form. Namely, assume (Xj, Yt)t>o 
is a one-parameter Gaussian semigroup which is not of the simple form 
flTU]) . Suppose that A, 5, H are as in Prop. |H We can clearly make an 
arbitrarily small change in the matrices A and H (hence getting new 
matrices A' and H') in a way that the set of eigenvalues of {A — H')a 
do not contain pairs of the type ±A. If necessary, we also make a 
small change to 5, obtaining S', to guarantee the condition iA' + B' > 
0. Hence, there exists a one-parameter semigroup (X/, Y^)t>o of the 
simple form which has generating matrices A',B', H' arbitrary close to 
A,B,H. 
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Example 3. Not all one-parameter Gaussian semigroups are of the 
simple form flTUl) . For instance, suppose n = 1 and choose A = 
together with 



H 



1 

1 



In this case 



A = {H- A)a 



-1 
1 



The generated matrix Xt = e^*^ corresponds to squeezing. 

With any y, the matrix Ay + yA is diagonal. Therefore, the con- 
dition f|T2|) cannot be satisfied whenever B (and hence B) is a non- 
diagonal matrix. On the other hand, since ^4 = any positive matrix 
is a possible choice for B. 

5.3. Bounded evolutions. We say that a one-parameter semigroup 
[Xt, Yt)t>o has bounded noise term if there is a constant c such that 
\\Yt\\ < c~for all t E M+. 

Proposition 6. Let {Xt, Yt)t>o be a one-parameter semigroup of Gaussian 
channels. The following conditions are equivalent: 

(i) There exists a positive matrix y such that 

Yt = y-XtyxJ. (13) 

(ii) [Xt, Yt)t>Q has bounded noise term. 

Proof. Suppose that (i) holds. For every t G M"*", we then have 

^<Yt = y-XtyxJ <y. (i4) 

Suppose that (ii) holds. Let /i be an invariant mean of the semigroup 
M'^ (see e.g. [I^). By the assumption, each matrix entry t i— > \Yi\ij is 
a bounded continuous function. We define the matrix y by 

As /i is invariant, an application to the second semigroup condition 
Yt+s = Yt + XtYgXf gives the formula (IT5]) . Moreover, 3^ > since 
Ft > for each t G M+ and 

{v\yv) = ii{{v\Y.v))>Q. 

□ 

Let {Xt, Yt)t>Q be a one-parameter semigroup and suppose there is a 
y satisfying (fTOj) . The matrix 3^ is a valid covariance matrix (e.g., of a 
Gaussian state) iff 

y-ia>Q. (15) 
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Thus, in this situation the one-parameter semigroup {Xt, Yt)t>o has an 
invariant Gaussian state. 

This is the case if the semigroup (XJ^^q is strictly contractive: 
< 1 for t > 0. Then we have limt^ooXt = so that ( |TT|) be- 
comes ffT5]) . 

Example 4. An attenuation channel is of the form X = y/rj 1, Y = 
(1 — r/) 1 for some t] G (0, 1). Attenuation channels form a bounded 
one-parameter semigroup. Namely, 

Xt = e-'l, Yt = (1 - e"^')l . 

For this one-parameter semigroup we have y = 1. Therefore, the 
vacuum state (with covariance matrix F = 1 and first moments d = 0) 
is an invariant state for the one-parameter semigroup of attenuation 
channels. 

We notice that a one-parameter semigroup {Xt, Yt)t>o may have 
bounded noise term without having an invariant Gaussian state. For a 
simple example suppose that {Xt, Yt)t>o is a one-parameter semigroup 
of reversible channels. Then Xt = e~^^'^ for some symmetric matrix 
H and Yt = 0. If Xt is not orthogonal (in which case 3^ oc 1) we have 
3^ = 0, which is clearly not a valid covariance matrix. 

6. Infinitesimal divisible channels 

We call a Gaussian channel parameterized by {X, Y) infinitesimal 
divisible if either 

(a) for every e > there exists a finite set of Gaussian channels 
{Xi, Yi) such that: 

(i) \\{X„Y^-{l,0)\\<e 

(ii) U^iX^,Yi) = {X,Y) 
or 

(b) the channel can be approximated arbitrarily well with (a)-type 
of channels. 

We note that in the classical case [5] this is what is called limit of a 
null triangular array, and it is proved there that such limits are precisely 
solutions of time-dependent Kolmogorov equations. Also in the quan- 
tum case solutions of time-dependent Markovian master equations are 
clearly infinitesimal divisible. For finite dimensional quantum systems 
the close relation between the two sets has been studied in [Ij. 

The concatenation of two infinitesimal divisible channels is clearly 
infinitesimal divisible. The infinitesimal divisible channels thus form a 
subsemigroup of Q. 
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It is clear that if {X, Y) is an element of a one-parameter semigroup 
of Gaussian channels, then it is infinitesimal divisible. The converse is, 
however, not true. Namely, suppose that S" is a symplectic matrix such 
that S ^ exp(sp(2n)). Then the channel (S*, 0) is not an element of a 
one-parameter Gaussian semigroup, but by Prop. [3]we have symplectic 
matrices Si, S2 such that S = S1S2 and Si, S2 G exp(sp(2n)). Since the 
channels (5*1,0) and (5*2,0) are infinitesimal divisible, so is also (S*, 0). 

The following is a simple necessary condition for a channel to be 
infinitesimal divisible. 

Proposition 7. If a Gaussian channel {X, Y) is infinitesimal divisible, 
then detX > 0. 

Proof. The claim follows from the continuity and multiplicativity of 
the determinant. □ 

Example 5. The Gaussian channel describing a phase conjugating 
mirror (with minimal noise) corresponds to the choices 



If n is odd, then det(X) < and therefore {X,Y) is not infinitesimal 
divisible. 

In the spirit of Corollary [T] we can formulate a converse of PropITl 

Proposition 8. If X & M2„(]R) satisfies det(X) > 0, then there exists 
a F e M2n(M) such that {X,Y) is an infinitesimal divisible Gaussian 
channel. 

Proof. Consider the real Jordan decomposition X = MJM~^ and 
group the Jordan blocks in J such that J = ® Jr with J„ be- 
ing the collection of all Jordan blocks with negative eigenvalues and 
Jr containing all the others. By defining Xi = M(— 1) © 1M~^ and 
X2 = M{-J_) © JrM-^ we get X = X1X2. Since det(X) > 0, the 
multiplicity of the eigenvalue —1 of Xi is even. By Proposition 1, the 
Xj's now have real logarithms which implies by Corollary [T] that there 
are l^'s such that {Xi, Yi) are elements of a one-parameter semigroup 
of Gaussian channels. Consequently, {X, Y) is infinitesimal divisible if 



We also have the following simple observation. 

Proposition 9. Let {X, Y) be infinitesimal divisible Gaussian chan- 
nel. Then every Gaussian channel {X, Y) with Y >Y is infinitesimal 
divisible. 




we choose Y = Yi + X1Y2XI. 



□ 
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Proof. We can split the additional noise Y — Y into arbitrarily small 
pieces, nr=i(l, - Y)) = {1,Y-Y). □ 

7. Divisible channels 

We call a Gaussian channel (X, Y) divisible if there exist two non- 
reversible Gaussian channels {Xi,Yi) and (X2,l2) such that 

{X,,Y,)-{X2,Y2) = {X,Y). 

In this section we show that actually all non-reversible Gaussian chan- 
nels are divisible — in spite of the existence of Gaussian channels which 
are not infinitesimal divisible. This is in sharp contrast to the finite 
dimensional case, where one has also indivisible channels [1]. 

The main tool which we use to prove the result is a surjective map- 
ping from Q onto the cone of positive matrices in M2„(C). Thus we 
define the mapping p in the following way: 

p : M2„(M) X M2„(M) M2„(C) , (X, Y) ^ t{XaX^ -a)+Y. 

A comparison of this definition with the cp-condition ([1]) shows that 
p{X, F) > iff (X, Y) e g. Moreover, p{X, F) = iff X G Sp{2n) and 
Y = 0. Thus, the kernel of p is exactly the set of reversible elements 
ofg. 

An essential property of p is the fact each positive matrix P is an 
image of a Gaussian channel {X,Y). This property is proved in the 
following two lemmas. This first lemma is a standard result in linear 
algebra, but we give a proof for the reader's convenience. 

Lemma 1. Suppose M G M2„(M) is anti-symmetric. Then it can be 
written in the form M = NaN'^ for some N G M2„( 



Proof. Since M is normal, there is an orthogonal real matrix R such 
that PTMR is a block diagonal matrix where each block is a 1-by-l 
matrix or a 2-by-2 matrix of the form 



a b 
—b a 



a,beR,b>0. (16) 



As M is anti-symmetric every 1-by-l block has to be and every 2-by-2 
block of the form flTBl) has a = 0. For this kind of 2-by-2 matrix we 
can write 

b\ f Vb \ f 1\ f Vb 



-b J Vov^A-io/Vov^, 

□ 
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Lemma 2. Let P e M2n{C) and P > 0. There ts {X,Y) G G such 
that p{X, Y) = P. 

Proof. Let us write P = Pr + iPi, where Pr,Pi G M2„(]R). Since 
P > 0, Pfj is symmetric positive and Pj is anti-symmetric. Also the 
matrix Pi + a is anti-symmetric, and by Lemma[T]it can thus be written 
in the form 

Pj + a = XaX^ , 

where X G M2„(M). Choosing Y = Pr the claim is proven. □ 

Each positive matrix P represents an equivalence class of Gaussian 
channels rather than a single channel. Namely, for two Gaussian chan- 
nels (Xi,F2) and (X2,F2) we have 

Yi) = p(X2, Y2) iff Fi = Y2 and X^crXf = X2(rX^ . 

If we concatenate two Gaussian channels corresponding to positive ma- 
trices Pi and P2 we obtain a Gaussian channel in the equivalence class 

P = Pi+XiP2Xf. (17) 

With this preparation we are now ready to prove that every non- 
reversible Gaussian channel is divisible. 

Proposition 10. Let {X, Y) & Q he non-reversible. There exist non- 
reversible {Xi,Yi), {X2,Y2) G g such that {X,Y) = (Xi,Fi) ■ {X2,Y2). 

Proof. Let us first consider the case detX = 0. Choose Xi = X, 
Yi = Y , X2 = 1 and I2 the projector onto the kernel of X. In this 
way we can write the channel (X, Y) as a concatenation of two non- 
reversible channels. 

Let us then suppose detX 7^ 0. Let P = p{X,Y). We denote 
Pi = eP, and fix a pair (Xi, Yi) such that p{Xi, Yi) = Pi. It is possible 
to choose < e < 1 such that det Xi 7^ 0. Indeed, by its definition Xi 
satisfies 

XiaXl = eXaX^ + (1 - e)a , (18) 

hence we get 

(detXi)2 = det(XiaXf) = det{eXaX^ + (1 - e)a) . 

For e = the right hand side is 1, and from the continuity of the 
determinant follows that det Xi 7^ for some < e < 1 . 
We then define 

P2:=(i-£)xr^pxr^. 

It follows that 

XiPaXf + Pi = P . 
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This shows that in the equivalence class of Gaussian channels repre- 
sented by P, there is at least one channel {X, Y) which can be di- 
vided non-trivially. So it remains to prove that this holds then for all 
Gaussian channels in the equivalence class. 

As we have noticed earlier, {X, Y) and {X, Y) are in the same equiv- 
alence class iff XaX'^ = XaX^. Thus, S := X-^X G Sp{2n) and 
X = XS^^. Therefore, a decomposition for {X,Y) also leads to a 
decomposition of {X,Y). 

Since Yi = sY ^ 0, the channel (Xi,Yi) is non-reversible. Also 
the channel (X2,F2) has to be non-reversible a.sYi ^Y implies that 
F2 ^ 0. □ 

Note that [Xi^Yi) can be chosen arbitrary close to the ideal chan- 
nel (1,0). That is, form an arbitrary irreversible Gaussian channel we 
can 'chop off' an infinitesimal (irreversible) piece so that the remaining 
part is still a valid Gaussian quantum channel. The possibility of iter- 
ating this procedure, i.e., chopping off an infinitesimal pieces from the 
remainder and so further, might suggest that every Gaussian channel 
is infinitesimal divisible. This intuition, however, fails since the re- 
maining channel (after having chopped off a piece) is, in the Gaussian 
context, not necessarily closer to the ideal channel (as it would be in 
the finite dimensional context). In fact, for Gaussian channel which 
are not infinitesimal divisible this procedure can bring us further and 
further away from the identity. This might be seen as a signature of 
the non-compactness of the set of Gaussian channels (as opposed to 
the compactness in the finite-dimensional context). 



8. Idempotent channels 

A Gaussian channel (X, Y) is idempotent if 

(X,F)-(X,y) = (X,F). 

This leads to the requirements X"^ = X and XYX^ = 0. Since Y is 
positive, the second condition can be written in the form {Xy/Y){XVYf = 
0, which is equivalent to XY = 0. Therefore, we conclude that {X, Y) 
is idempotent iff 

X^ = X, XY = 0. (19) 

In physical terms, idempotency means that a repeated use of the chan- 
nel does not change the system any further. 

Suppose that X is invertible. Then the conditions (fT9!) imply that 
{X,Y) = (1,0), which is just the identity channel. However, there are 
also other idempotent channels as illustrated in the following example. 
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Example 6. Let X and Y be diagonal matrices of the form 

2k 2n-2k 2k 2n-2k 

X = diag(l, . . . , 1,0, . . . ,0) , Y = diag(0, . . .,0,yi,yi, . . . ^yn-k,yn~k) , 

where yj > 1; j = 1, . . . ,n — k. Then the pair (X, Y) clearly satisfies 
f ll9p . The cp-condition ([T]) breaks into conditions for 2 x 2-matrices, 

iai > iai , 2/j 1 > ^o"i , 

which obviously hold. The channel (X, Y) corresponds to a transfor- 
mation where we do nothing for the first 2k modes but for the rest 
2n — 2k modes we do a state preparation (see Example [1]). 

Let S G Sp{2n). If we concatenate a channel {X,Y) with the re- 
versible channels corresponding to 5* and 5*"^, we get 

{S, 0) ■ (X, Y) ■ {S'\ 0) = {SXS-\ SYS^) . (20) 

It is easy to verify that if {X, Y) is idempotent, then also {SXS~^, SYS'^) 
is idempotent. Therefore, Example [HI generates a full class of idempo- 
tent channels. In the following we show that actually all idempotent 
channels are of that form. 

Proposition 11. A Gaussian channel {X,Y) is idempotent iff there 
is a symplectic matrix S such that 

2k 2n-2k 

SXS-^ = diag(l7'^, 0^^^~0) , (21) 

^ 2n-2fc 

SYS'^ = diag(0, . . . , 0, yi, yi, . . . , yn-k, y-a-k) , (22) 
where yj > 1; j = 1, . . . ,n — k. 

Proof. Let us first show that X^ is a symplectic projection, i.e. the 
symplectic space (M^", a) is a direct sum of two subspaces Vi, V2, mu- 
tually orthogonal with respect to the symplectic form. That is every 
vector V is uniquely decomposed as vi + V2 with vi G Vi,V2 G V2, and 
vjav2 = 0. 

Indeed, put Vi = X'^v,V2 = {1—X'^)v, then v = V1 + V2 and yt^i = 0. 
By condition ([T]) we get 

V2YV2 = {vi + iv2yY{vi + iv2) > 2v'[av2 

for all vi E Vi,V2 E V2, hence v'[av2 = 0. 
From ([T]) we also obtain 

Y = {1- X)Y{1 - Xf > 2(1 - X)a{l - Xf . 
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Applying symplectic diagonalization of the positive symmetric matrix 
(1 —X)Y{1 — X)-^ in V2, one can always find a symplectic matrix 5* in 
(i?2",cr) satisfying 
Since XY = 0, we also have SXS~^SYS^ = and ([21]) follows. □ 

9. Gauge- covARiANT channels 

We say that a Gaussian channel {X, Y) is gauge- covariant if [X, a] = 
[Y, a] = 0. We denote by Qa the set of all gauge-covariant Gaussian 
channels. It is clearly a subsemigroup of Q. Physically, gauge-covariant 
channels arise for instance from a number conserving (i.e., passive) 
coupling to an environment (cf. |17j). 

Let us rearrange the matrix a such that 

1„ 
-1„ 



cr 



then matrices M G M2„(C) commuting with a are those of the form 

where A,B ^ M„(C). The maps M A ± iB are easily seen to 
be *— homomorphisms of M2„,(C) onto M„(C), hence M > implies 
A±iB > 0. Note that A* = A,B* = -B in this case. Let us 
show that conversely, M > ii A ± iB > 0. Let A ± iB > 0, then 
A > and B* = —B. Assume first that A is nondegenerate, then 
l„±zA-i/2^y4-i/2 > 0^ which implies ln + A-^^^BA^^BA~^^^ > and 
hence A + BA~^B > which implies M > 0. The case of degenerate 
A is obtained by approximation. Thus we have proved 

Lemma 3. M > iff A ± iB > 0. 

For a matrix M of the form (!23|) we denote M = A + iB. Then 
a = iln and the cp-condition ([1]) for gauge-covariant channels takes 
the form 

Y >±(ln- XX*] . (24) 



Let X* = UK be the polar decomposition of the matrix X*, where U 

is unitary and K = Vxx* is positive. Then the channel {X, Y) is a 
concatenation of the reversible channel {U,0) and the channel {K,Y). 
For the last channel the condition (l24l) takes the form 

f > ± (l„ - . (25) 
There are several basic cases depending on the properties of K: 
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(i) K = 0. The channel is idempotent for any Y > 1„. (State 
preparation with the covariance matrix Y). 

(ii) < K < In, where by strict inequahty we mean that the eigen- 
values kj of K satisfy < kj < 1. The channel is a member of 
a one-parameter semigroup of Gaussian channels with bounded 
noise term having an invariant state. The semigroup is defined 
via the relations Xt = K*, Yt = y — K^yK*, where y is the 
unique solution of the equation 

Y = y-kyk 

corresponding to the covariance matrix of the invariant state. 
Indeed, the last equation written in the basis of eigenvectors of 
K has unique solution 

y = (uij) , Uij = (1 - kikj)-^ yij, Y = {jjij) . (26) 

Let us prove that 3^ > 1„ and hence it corresponds to the 
covariance matrix of a Gaussian state. The matrix with the 
elements 

roc 

(1 - kikjY^ = / expt {hkj - 1) dt (27) 
Jo 

is positive since the matrix {kikj — 1) is conditionally positive 
(see Thm. 6.3.6. in [IS]). By using the condition 1" > 1„ — 
and a lemma concerning Schur products of positive matrices 
(Thm. 5.2.1. in [16]), we have 

y={{l- k,k.j)-^ y,,) > ((1 - k,k,)-' (1 - k^) 6,,) = 1„. 

(iii) K = In- Then Y > and the channel is member of one- 
parameter semigroup of Gaussian channels {l,tY) with un- 
bounded noise term. In the case Y = this is identity channel. 

(iv) K > In, that is kj > 1. The channel is a member of one- 
parameter semigroup of Gaussian channels with y defined as in 
fl2Bl) . However in this case instead of fl^7|) we must use 

POO 

(1 - kikj)'^ = - expt (1 - kikj) dt, (28) 
Jo 

therefore the matrix with the elements (128|) is negative implying 
3^ ^ —In- Thus the semigroup has unbounded noise term and 
there is no invariant state. 

In general, one can decompose K into direct orthogonal sum of the 
matrices satisfying the conditions (i)-(iv). In case Y commutes with 
k one can further decompose the channel {K, Y) into corresponding 
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channels. In particular, in the case n = 1 any gauge-invariant channel 
reduces to one of the cases (i)-(iv). 



10. Conclusions and open questions 

We conclude with some open questions. First of all, we lack a charac- 
terization of infinitesimal divisible Gaussian channels. ProplH] provides 
a partial answer in terms of the determinant det(X). A similar prop- 
erty, however, turned out to be necessary but not sufficient in the case 
of finite dimensional quantum channels |1] . Moreover, we left open the 
question which type of dynamical equations (e.g., time-dependent Mar- 
kovian master equations) leads to solutions which coincide with the set 
of infinitesimal divisible channels. 

For one-parameter semigroups our picture is more complete. Yet, 
there is no efficiently decidable criterion which enables us to say whether 
or not a given Gaussian channel is an element of such a one-parameter 
semigroup. The simple form (fTOj) suggests to follow the lines of [5l [6] 
where an integer semi-definite program provided a solution for finite 
dimensional quantum channels. However, boundary cases (e.g., chan- 
nels not admitting a simple form representation) will have to be treated 
with care. 

Other questions arise when we slightly change the rules of the game. 
In the reversible case we saw that while a transformation might not be 
an element of a one-parameter semigroup within the Gaussian world, 
it can become one if we drop the restriction to the Gaussian world. So 
how does the general picture change if we allow for decompositions into 
arbitrary channels? 

In a similar vein we may allow for tensor powers and thereby in- 
vestigate the robustness of all the discussed properties w.r.t. taking 
several copies of a quantum channel. This might be interesting beyond 
Gaussian channels (e.g., for qubit maps) as well. In the Gaussian case 
we can easily find examples showing that things can change: take a 
reversible Gaussian channel with 5* ^ exp(sp(2?7,)), then S (BS happens 
to have a Hamiltonian matrix generator. 

Finally, it is desirable to relate semigroup properties of a quantum 
channel to other properties such as their capacities or to properties of 
quantum spin chains to which the channels can be assigned to via the 
finitely correlated state construction [18], [19]. 
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